Where is the Luttinger liquid in one dimensional semiconductor quantum wire 

structures? 
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We present the theoretical basis for analyzing resonant Raman scattering experiments in one- 
dimensional systems described by the Luttinger liquid fixed point. We make experimentally testable 
predictions for distinguishing Luttinger liquids from Fermi liquid and argue that presently available 
quantum wire systems are not in the regime where Luttinger liquid effects are important. 
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It is theoretically well-established [1-3] that a one di- 
mensional interacting electron system (IDES), is not a 
Fermi liquid (FL). Unlike a Fermi liquid, the interact- 
ing IDES has neither sharp fermionic quasiparticle exci- 
tations nor a discontinuity in the electron momentum 
distribution function. The elementary excitations are 
charge e, spin zero bosons and spin 1/2 charge 'semions' 
(fractional statistics objects), and the fermion is a com- 
posite of these. Interacting IDES have been generically 
termed Luttinger Liquids (LL) [2] and have been the sub- 
ject of extensive theoretical study over the last 40 years 
and particularly over the last decade. Despite the intense 
theoretical interest, there have been few convincing ex- 
perimental demonstrations of the predicted LL behavior 
in real IDES. The power-law density of states observed 
in tunneling into edges of quantized hall systems [4] have 
been interpreted in terms of the theoretically expected 
'chiral Luttinger liquid' behavior of edge states. The ori- 
gin of the differences between the observed and expected 
exponents is presently an area of active inquiry. Pho- 
toemission experiments on Mott insulating oxides have 
been interpreted in terms of the 'holon' and 'spinon' ex- 
citations of a charged Luttinger liquid [5]. 

A IDES which is of particular interest both for funda- 
mental physics and for technology is the system formed in 
GaAs-based semiconductor quantum wire (QWR) struc- 
tures. Modern materials growth and fabrication tech- 
niques have produced nearly ideal IDES in which the 
electron may move freely only along the length of the 
wire. The transverse motion is quantized with the quan- 
tum ID subbands separated by several meV. It is pos- 
sible to have low enough carrier densities so that at low 
temperatures only the lowest ID subband is occupied by 
electrons. Such GaAs QWR based IDES should be ideal 
systems for the study of interacting electrons in one di- 
mension because they are free from complications aris- 
ing from band structure, lattice effects, and crossovers to 
three dimensional behavior which often make interpre- 
tations of experimental data difficult in more traditional 
IDES based on organic compounds. 

It is surprising, therefore, that no definitive LL be- 



havior has been reported in GaAs QWR systems, and 
in fact the ID Fermi gas/liquid model seems to "work" 
operationally very well in describing and explaining the 
observed IDES experimental properties in GaAs QWR 
[6,7]. Part of the reason for the apparent absence of the 
expected LL behavior is undoubtedly the fact that in 
weakly interacting IDES, at finite temperatures and in 
the presence of impurity scattering, the actual quantita- 
tive difference between a LL and a FL is not large [7], 
although the qualitative conceptual difference between 
the two is huge. A more fundamental issue is that the 
differences between a Luttinger Liquid and a Fermi liq- 
uid are most obvious in the one-electron spectrum, while 
the experimental probes which may most conveniently 
be applied to the QWR structures produce 'particle-hole 
pairs'. The differences in particle-hole pair properties 
between Luttinger Liquid and Fermi Liquid systems is 
much less pronounced than are the differences in the one 
electron spectrum. This perhaps accounts for the fact 
that one of the most important probes of QWR struc- 
tures, resonant inelastic light scattering or Raman scat- 
tering spectroscopy (RRS) [6,8], has not yet observed any 
definitive indications of LL behavior in these systems. 

In RRS experiments, light is absorbed at one frequency 
and re-emitted at another, creating one or more particle- 
hole pairs. In the so-called polarized geometry with the 
incident and outgoing photons having the same polariza- 
tion (so that no spin is transferred to the QWR), RRS 
experiments in GaAs QWRs consistently [6,8] show two 
peaks which indeed look qualitatively very similar [9] to 
the spectra for the corresponding 2D and 3D systems. 
In these higher dimensional systems, the two peaks have 
a clear and generally accepted Fermi liquid interpreta- 
tion [9]. The higher energy peak is associated with the 
plasmon or charge density excitation (CDE), a collective 
density excitation of the electron gas, and the lower en- 
ergy spectral peak is associated with incoherent particle- 
hole pair excitations (SPE). In the QWR materials, the 
lower energy peak occurs at an (approximate) excitation 
energy of ~ qvp, where q is the excitation momentum 
and vf is the ID Fermi velocity obtained from the band 
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structure of the QWR. An interpretation of the lower 
peak as an SPE contribution seems therefore natural [9] . 
However, there is a strong theoretical objection to this 
interpretation: in a one dimensional system there is spin- 
charge separation: the only charge excitations live at the 
plasmon frequency, and cannot contribute to excitations 
at the SPE energy. The signal observed in this q, v range 
must be due to the chargeless spin excitations of the LL; 
in particular it is possible to combine two S — 1/2 exci- 
tations into a, S — object, creation of which is allowed 
by the Raman selection rules. 

Sassetti and Kramer (S-K) presented a qualitative the- 
ory of this effect [10]. They showed that although the 
leading contribution to the RRS matrix clement corre- 
sponds to coupling the light to the electron density op- 
erator, there is a sub-leading term (which becomes more 
important under resonance conditions) which may be in- 
terpreted as a coupling of light to the energy density 
fluctuations of the electrons in the QWR. The energy 
density fluctuations have a contribution from the spin 
excitations, which qualitatively explains the data, but 
the S-K theory did not calculate the spectral weights of 
the RRS peaks. Too close to resonance, the S-K the- 
ory breaks down. The S-K work also does not show how 
to distinguish a LL from a FL in the RRS experiment. 
The most important theoretical problem is that the S- 
K calculation is logically inconsistent, because it uses an 
expression for the RRS matrix element which is correct 
only if the conduction band is a FL not an LL. Thus S- 
K uses FL matrix elements but LL excitations. In our 
paper the correct LL matrix element is used, leading to 
expressions different from those derived by S-K. 

In this paper we present an essentially complete treat- 
ment of RRS in a one dimensional electron gas. We ob- 
tain a precise expression for the energy transferred to 
the QWR in a RRS experiment, valid at all values of the 
difference of the energy from resonance, and evaluate it 
quantitatively in several experimentally relevant limits. 
We show which features of the data contain information 
about the LL exponents, obtain expressions for the rel- 
ative amplitudes of the SPE and CDE peaks, determine 
lincshapes and discuss qualitatively the crossover from 
LL to FL behavior. 

Resonant Raman scattering is a two-photon process in 
which a photon is absorbed, transferring an electron from 
the valence (V) band to the conduction (c) band and a 
photon is emitted, transferring an electron from the con- 
duction band back to the valence band. We assume that 
the valence band is initially filled, and assume there is 
no excitonic interaction between conduction and valence 
band states. The excited valence hole is then described 
by a single-particle Hamiltonian, which we write as Hy, 
while the conduction band is described by some interact- 
ing Hamiltonian which we denote Hll- We denote the 
photon absorption and emission by respectively. The 
RRS process is described by the following Hamiltonian: 

H = H V + H LL + P 1 +P 2 (1) 



where the photon- in (Pi) and photon-out (P2) terms are 
Pi = e-^+^X £4 +9 /2, s (*K*(*) +h.c. (2) 

p.s 

p 2 = e i ( n-„/2)t £ vl s (t)c P - q/ 2At) + h-c. (3) 

with c and v the annihilation operators for electrons in 
conduction and valence band states respectively. Note 
that the operator vj, CT creates an eigenstate of Hy with 
energy E% while the cj, i<7 operators does not create eigen- 
states of Hll- The absorbed (emitted) photon energy 
and momentum are set VL ± v/2 and ±q/2 respectively. 

We now use the standard methods of time-dependent 
perturbation theory to calculate the amplitude, a n (to), 
for the system at some time to to be in some excited 
state |n) of QWR, but with no holes in the valence band. 
We assume the system is in its ground state at t = 0. Our 
neglect of any excitonic interaction between conduction 
and valence band simplifies the calculation and we obtain 

a n (to) = T V / dRer* 11 [° dTe l » T (n\d rs (R,T)\0) 
L r,s J J o 

(4) 

with 

d rs (R,T) = J dx J dt<f>(x,t) x 

As(R + x/2, T + t/2)4 s (R - x/2, T - t/2), (5) 
where r and s are band and spin indices (pml), and 

(P(x,t)=e int J2 ei{E ^~ px) - ( 6 ) 
p 

Eqs. (4) and (5) are our fundamental new results: they 
show that the RRS process acts to create a particle-hole 
pair at a spatial separation x and temporal separation t. 
These arc determined by the average photon frequency £1 
and the valence-band properties encoded in E p . Further, 
if interactions are present in the conduction band, the 
states created by ^ and by tp are not eigenstates of Hll 
and therefore the matrix element is itself modified by 
interactions. 

We note that Eqs. (4) and (5) maybe substantially 
simplified in the limit of greatest physical interest. We 
linearize the valence band energy about the conduction 
band Fermi momentum, writing Ep = — A — Vp (rp — pp) 
for branch r and define ujr = ft — A as the photon fre- 
quency with respect to the resonance energy, A. The 
p— integral gives 8{x + Vpt). Finally we write the conduc- 
tion band operators in terms of the bosons which create 
eigenstates of Hll, and normal-order in the boson basis, 
obtaining 

O rs {R,T) = L [ dte lbjRt G c rs (-rv^t,t) 
Jo 

. p i®rs,i,(R,-rv^t;T,t) .. i$ re ,„(R,-rv%t;T,t) . 
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FIG. 1. Calculated polarized RRS spectra for various reso- 
nance condition, lor. One- and two-boson contributions have 
been plotted separately in order to show their relative con- 
tributions (see text). A finite broadening 7 has been used 
to depict the results. Note that the overall spectral weights 
decreases dramatically off-resonance, as indicated by the in- 
dividual scale factors on right side of each plot. 
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FIG. 2. Spectral weights for the low energy (solid curve) 
and the high energy (dashed) RRS peaks in the LL theory 
plotted A function of the Luttinger exponent a: lor — 
(main); 0.1 (inset). When \lor\ larger than 0.1 the low en- 
ergy ("SPE") weights are always much smaller than the high 
energy (CDE) weights over the whole range of a. 
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p>0 



{ - sinh 6 P sm[p(rx + v p t) /2] [bl rp e lp(rR+v ? T) + h.c] 
+ cosh6 p sm[p(rx - v p t)/2}[bt p e- ip(rR -^ T) + h.c.]} 



(8) 



* ra>a {R,x;T,t) = 2sy e- ap l\ ^- x 
P >o V P L 

{sm[p(rx - v c F t)/2][al p e-^ rR - vC ^ + h.c.]} . (9) 
Here b + and er + create charge and spin excitations respec- 



tively and v p 



is the plasmon velocity, where 



the exponent e~ 29p = y/\ + 2g/nv F is defined for the 
short-ranged interaction, g. G c is the exact conduction 
band Green's function at spatial separation —rv F t, and 
temporal separation t. We have assumed that the in- 
teractions are negligible in the spin sector and therefore 
the spin excitation velocity is just the Fermi velocity. 
As long as v F , the valence band velocity at the conduc- 
tion band pf is different from the spin and charge veloc- 
ities of Luttinger liquid, G c is a decaying function of t. 
In the noninteracting case, G c ~ 1/t; interaction leads 
to a faster decay: G c ~ l/t 1+a with the LL exponent 
a = sinh 2 6 p > (not the same one as we use in Eqs. (8- 
9) for infinitely small convergent factor) for short-ranged 
interactions; G c decays faster with the physically rele- 
vant long-ranged interactions. This faster decay of G c 
is the mathematical expression of the renormalization of 
the RRS vertex by the interactions, which produce the 
Luttinger liquid behavior. As we will now show, it has 
important consequences for various aspects of the RRS 
spectra; and in particular for the dependence of the CDE 
and SPE energies on the difference of the average photon 
energy from the resonance. 

We defer to a subsequent paper a full evaluation of 
the RRS correlation function, which is computationally 
demanding and not very illuminating, and present here 
the results of expanding Eq. (7) in terms of boson oper- 
ators. The essential point is that if the combination of 
e tuJat G c (—rv F t,t) decays rapidly as t increases (large lor 
as off- resonance or large a as strong interaction), then 
the t— integral is dominated by small times and an ex- 
pression in power of bosons is rapidly convergent. We 
will show below that the first order term, one-boson re- 
sult, gives the main contribution to CDE spectrum and 
dominates the whole spectrum as off-resonance and the 
second order term, two-boson (spinon) result, gives the 
peak at "SPE" energy as near resonance, but it still has 
relatively small weights as compared to the first order 
CDE. 

Expanding the exponentials, keeping only the one- 
boson term and integrating explicitly, gives the one- 
boson transition rate as a delta function at v = qv p with 
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the spectral weight (a < 1) 



Wx 



2LT 2 (-a) 



Eq 



( U R 



V E 



(10) 



qVp/2, neglecting Vp for simplicity. Eq is 



where u> 

the energy scale depending on the interaction range and 
roughly of the order of Fermi energy, E C F . For uj q <C \wr\, 
W x oc M 2a ~ 2 , while for w R = 0, W x « sin 2 (7ra/2). 
Thus LL effects enter the CDE portion (one-boson) of 
the spectrum in two ways (for short-ranged interaction) : 
first, far from resonance, it changes the frequency de- 
pendence of spectral weights from w^ 2 , the noninter- 
acting result, to Wr a (note that all other higher or- 
der bosonic contribution decays much faster, this con- 
firms the validity of the bosonic expansion we mentioned 
above), and secondly as on resonance (ojr = 0) it changes 
the value to be nonzero due to finite interaction strength. 

To second order, two new effects appear. In the density 
spectrum, branch mixing process appear. These lead to a 
continuum absorption beginning at the CDE threshold, 
lu = qv p . In addition, an S — combination of spin 
excitations may be excited via the two spinon, (aa, aa) 
(note that there is no first order contribution in spin 
channel due to the selection rule of RRS in the polarized 
spectroscopy), and gives the so-called "SPE" mode at 
V = qvp. 

In Fig. 1, we show the spectrum from one and two 
bosons in different resonance energy. One can find that 
(i) the overall spectral weights decays very fast off reso- 
nance, and (ii) the "SPE" peak is generated at u> ~ 0.2Ep 
by the two-boson contribution near resonance. But as 
compared with the CDE peak at plasmon energy (about 
0.57 Ep), the "SPE" peak is still very small compared 
with CDE. This striking result arises from the fact that 
the contribution of one spin-boson in the first order is 
forbidden by the specific selection rule of polarization 
in depolarized RRS spectroscopy, (iii) At higher energy 
side above CDE peak, there is some continuum structure 
which is not shown in the range of Fig. 1. This contin- 
uum is from the interaction between different branches 
of charge bosons due to finite gi interaction. We are not 
interested in their structure because it goes to zero near 
the plasmon energy and their higher energy behavior is 
off the experimentally measurable region, and become 
unphysical due to the failure of the linear dispersion as- 
sumption, (iv) When including three or higher order bo- 
son contribution (not shown in this paper), we will see 
the mixture of charge boson and spin boson in a form like 
{crap, po~o~), which will come into the energy between qvp 
and qv p , plasmon energy, as a continuum structure. A 
detailed analysis shows that these higher order contribu- 
tion is relatively small and no special structure compared 
to the first two order result we present here. While Fig. 
1 is for a specific value of a (— 0.3) we show in Fig. 2 
the calculated charge boson and spin boson RRS spectral 
weights at resonance and away from resonance. In gen- 
eral, the LL theory predicts much smaller spectral weight 



for the lower energy " SPE" mode than the FL theory [9] 
at resonance. This is particularly true since our best es- 
timate for the Luttinger exponent of the experimental 
system [6] (obtained from the CDE energy dispersion) is 
a ~ 0.4. 

As compared with the experimental result, which 
shows possible comparable spectral weight of " SPE" with 
CDE [6] , we find that the LL theory result induced by res- 
onance effect does not explain the experimental results 
quantitatively, even though we could recover the SPE 
peak through the coupling of two spinon in LL, not as 
the SPE in FL theory. This inconsistency cannot be re- 
solved even evaluating the full bosonic contribution with- 
out expansion as what we have in this paper. Our future 
work shows that the spectral weight of the " SPE" peak 
enhanced by spinon coupling in polarized RRS spectra is 
always relatively small compared to that of CDE. There- 
fore, unlike the conclusion of previous work [10] based 
on the incorrect matrix element, we claim that the whole 
problem cannot be simply understood by the correct LL 
theory. We believe that the existing experimental results 
[6] are in the high energy crossover regime where in fact 
a FL description maybe more appropriate for the RRS 
data than the LL description which is an asymptotic low 
energy description. This explains the spectacular quan- 
titative success of the FL RRS theory developed in ref. 
[9]. 

In conclusion, we provide the correct LL theory for 
the RRS spectra calculation, and obtain some meaning- 
ful and interesting results to study the possible origin 
of LL features in the RRS spectra of ID QWR systems. 
We also develop an useful bosonic expansion method to 
study the two-particle correlation function. Finally, we 
find that the LL theory cannot quantitatively explain the 
experimental data most likely because the RRS experi- 
ments are not in the asymptotic low energy LL regime. 
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